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Abstract 



A one-parameter family of variational problems is introduced that interpolates between tor- 
sional rigidity and the first Dirichlet eigenvalue of the Laplacian. The associated partial 
differential equation is derived, which is shown to have positive solutions in many cases. 
Results are obtained regarding extremal domains and regarding variations of the domain or 
the parameter. 

1 Introduction 

Various geometric and physical constants can be associated with a bounded domain in the plane: 
perimeter, area, diameter, transfinite diameter or capacity, torsional rigidity and principal fre- 
quency. We focus on the latter two physical quantities. The torsional rigidity P(D) of a bounded 
domain D in the plane is given by 



Each has a physical significance. In case D is simply connected, the torsional rigidity is a measure 
of the torque required per unit length per unit angle of twist when a beam with cross section D 
undergoes torsion. The principal frequency is a measure of the lowest note that a drum of shape D 
can produce. An extensive literature is devoted to each. The classic book |PSj by Polya and Szego 
contains a wealth of information on both, in particular the variational characterizations ([I]) and 
|2]) and the Euler-Lagrange partial differential equations (PDEs) listed below. The corresponding 
Euler-Lagrange equations for these two problems are 




(1) 



and the principal frequency A(D) is given by 




Au + 2 = 0, 




(3) 



for torsional rigidity and 

Au + Xu = 0, u\ QD = (4) 



for frequency. Classical theorems tell us that both these equations have positive solutions for a 
bounded domain D, provided the boundary 3D has some regularity. 

One can measure the fundamental frequency of a bounded domain in any dimension; both the 
variational formulation ([2| and the PDE Q remain the same. The torsion function, the solution 
to the problem ([3]), is of probabilistic significance in any dimension, being the expected exit time 
of Brownian motion from the domain D relative to Wiener measure. 

The purpose of this note is to introduce a one-parameter family of variational problems and 
corresponding PDEs that interpolate between the torsional rigidity and principal frequency prob- 
lems listed above, and go somewhat beyond them. In the remainder of this introduction, we briefly 
describe the interpolating problem and the relevant physical quantities, state some results, and 
outline the remainder of the paper. 

Let D C R n , n > 2, be a bounded domain with smooth boundary (i.e. one can locally write 
the boundary of D as the graph of a C°° function). 

Definition 1. For each p > 1, we define C p (D) by 

C p (D) = inf I = |V " |2 ^ p :ueLP(D)n W^ 2 (D),u^ o\ . (5) 



{J D uPdV) 



In particular, A/C\(D) is the torsional rigidity and 02(D) is the fundamental frequency. The 
scaling law for C p is 

C p (rD) = r n - 2 -^C p (D), r > 0. (6) 
In Section [2] we derive the corresponding Euler-Lagrange equation. 

Theorem 1. Let D C R" be a bounded domain with smooth boundary. Let p > 1. The critical 
points of the functional 



r l,2/ 



(f D v?dV) 2/p 

in L P (D) n Wq'*(D) satisfy the PDE 

Au + Au p - 1 =0, u\ dD =0 (7) 

for some constant A. 

After verifying that this is the correct Euler-Lagrange equation, we discuss solvability of the 
PDE, borrowing some classical results of Pohozaev [Pohj . and some basic properties of solutions. 

In contrast to the functional $ p , the differential equation ^ is not scale-invariant unless 
p = 2. If u solves Q for some constant A and k > is a constant, then v = ku satisfies 

Av + k^AvP- 1 =0. 

Conversely, given a solution u to Q and a constant a > 0, we see that 

v = ^— ^ P u solves Av + a/ -1 = 0. 

Thus, if p ^ 2, by rescaling we obtain solutions to the equation Au + Au^" 1 = for all A > 0. 
The variational problem, however, will select a particular Lagrange multiplier A. 

Lemma 2. Let C p (D) be given by ^ and let u be a positive solution to ([7| for some A > 0. 
Then 

/ r \ (P~ 2 )/P 

C P (D) = A ( J u p dVj . (8) 
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Notice that the right hand side of (|Sj) is invariant under scaling of the function u. We recover 
a classical formula for the torsion function when p = 1, namely P{D) = 2 J D udA where u is the 
solution of ([3]) . The case p = 2 returns the first Dirichlet eigenfunction of the Laplacian. 

Motivated by this lemma and torsional rigidity, we make the following definition. 



2n 



Definition 2. Let D c K." be a bounded smooth domain and either n — 2, or n > 3 andp < . 
Then the p-torsional rigidity is given by 



K p (D) = ~UuVdV 



(2-p)/p 

(9) 



where u > solves the boundary value problem |?j). 

Then Tl x {D) = P{D) and K 2 (D) = 4/A(D). We see immediately from the definition and ^ that 

C p {D)K p (D) = l, p>\. (10) 

Next, in Section |3j we examine the specific case where D is a round ball or an infinite slab 
{ — 1 < x n < 1} C K™. We prove some comparison results for C P (D) in Section|4j varying either 
the domain D or p. In particular, we prove the following theorem. 

Theorem 3. Let D C R™ be a bounded domain with volume V — Vol(D), and let 1 < p < q. 
Then 

V 2/p C p {D) > V 2/q C q {D). (11) 



In dimension 2, the case p — 1 and q = 2 of the inequality (11) becomes X(D) < 4A(D)/ P(D) 
which one may find on Page 91 of |PS] and relates the fundamental frequency of a domain to its 
area and its torsional rigidity. The scaling law ^ shows that 

V(rD) 2 / p C p (rD) = r n - 2 V{D) 2 / p C p (D), p>l, (12) 

which agrees with the classical scaling laws for the torsional rigidity and the principal frequency. 
In Section [5] we characterize extremal domains for C p under certain conditions. We prove an 
inequality of Faber-Krahn type. See (0 |K] for the classical Faber-Krahn inequality for principal 
frequency (the case p — 2) and see [P] or Appendix A of |PSj for Polya's proof of the Saint- Venant 
theorem that among all simply connected domains of given area the disk has the largest torsional 
rigidity. 

Theorem 4. Let p > 1, let D C M. n be a bounded domain with C°° boundary, and let B C W 1 
be the round ball, centered at the origin, with the same volume as D. Then C P (D) > C P (B). 
Moreover, equality only occurs if D = B almost everywhere. 

We also discuss convex domains of fixed inradius in this section. We conclude with a short list of 
open question in Section [6j 



2 The variational problem and its Euler-Lagrange equation 



We take p > 1 and a bounded domain DcR™ with a C°° boundary, and for u € L p (D)C\Wq 2 {D) 
not identically zero define the functional 

_ J D \Vu\ 2 dV \Wu\\ 2 L2{D) 

pi ) ~ (J DUPd vf p ~ MU D) ■ 

Our first task is to derive the Euler-Lagrange equation 



3 



Proof. Observe that $ p is scale-invariant; that is, if k > then = $ p (it). Thus, we can 

reformulate the condition that a is a critical point of $ p as a constrained critical point problem: 
find the critical points of J D \Vu\ 2 dV subject to the contraint J D u p dV = 1. Any constrained 
critical point must satisfy 



d_ 
de 



\V(u + ev)\ 2 dV = A 



e=0 J D 



(u + evfdV 



for all v £ L P (D) n Wj' 2 (D), where A is the Lagrange multiplier. Next recall that C^°(D) is dense 
in both L P (D) and W ' (D) (see, for example, Section 7.6 of |GT| ). so without loss of generality 
we can take u, v £ C^(D). Thus we can freely differentiate underneath the integral sign, and a 
quick computation shows that the equation above is equivalent to 

0= f -(Vu,Vv) +Au p - 1 vdV= [ v{Au + Au^dV. (13) 

JD JD 

Here we have absorbed a factor of 2 and a factor of p into the Lagrange multiplier A. If equation 
(131 is to hold for all compactly supported v in D, then u must satisfy the PDE 

Au + AuP- 1 = 

for some constant A as claimed. □ 

Remark 1 . This is a familiar differential equation, occurring (for instance ) in the study of scalar 
curvature under a conformal change of metric; see \LPf . 
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Remark 2. One can equally well study the functional 

(j D \Vu\ W) 
(f D v?dV) 2/p 

In this case, the Euler- Lagrange equation is 

= A q u + Au 13 - 1 = div(|Vu|«- 2 Vu) + AuP' 1 . 

This differential equation is either singular (for q < 2) or degenerate (for q > 2) at the critical 
points ofu. Thus, we do not expect to have as well-developed a theory attached to the more general 
variational problem. 

To examine the solvability of equation Q, we first recall the following classical theorem of 
Pohozaev |Poh] . 

Theorem. (Pohozaev) Let D C K n be a bounded domain with smooth boundary, and let f(u) be 
a Lipschitz function. If n = 2 and f satisfies the estimate 

\f(u)\ < A + B\u\e c ^\ a<2, 

then one can find eigen] 'unctions of the PDE 

Au + \f(u) = 0, u\ 9D = 0. 

If n > 3 and f satisfies the estimate 

n 4- 9 

\f(u)\ < A + B\u\ m , m< - 



n-2' 
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then again one can find eigenf unctions of the PDE 

Au + Xf(u)=0, u\ gD = 0. 
Conversely, if D is star-shaped with respect to the origin and u > 0, not identically zero, solves 

Au + u m = 0, u\ gD = Q 

then m < . 

n—2 

Applying Pohozaev's theorem, we immediately obtain the following corollary. 

Corollary 5. There is a positive solution to the Euler- Lagrange equation if either n = 2, or 
n > 3 and p < ^rg. On the other hand, if n > 3, p > ^3? , and D C W 1 is star-shaped, then 
equation does not have a positive solution. 

It is well-known that for the critical value p = a minimizing sequence for the function $ p 
will typically become unbounded. Thus it is difficult to determine whether one can realize the 
corresponding infimum C P (D) as & P (u) for some u > 0. One can find a treatment of this blow-up 
phenomenon, which reflects the loss of compactness in the Sobolev embedding theorem, in |Tr| 
and Section 4 of [LP] . 

We use the maximum principle to prove the following lemma. 

Lemma 6. Let D C K n be a smooth, bounded domain and A a constant. Then there is at most 
one positive solution to the boundary value problem 

Proof. Suppose u and v are distinct positive solutions to ([7]), and let D* be a connected component 
of {u ^ v}. Without loss of generality we can assume that A > and that u > v on D* . We can 
also assume dD* is smooth by Sard's theorem. Then u — v is a positive solution to the boundary 
value 

A(u -v) = AK- 1 - v*- 1 ) > 0, u\ dD , = v\ dD , , 
which contradicts the maximum principle. □ 

If D is a convex domain, we can in fact say more. A theorem of Korevaar (see Theorem 2.5 
of |Kor| and the remark immediately following it) implies 

Corollary 7. If D C K." is a bounded, smooth, strictly convex domain and u > solves the 
boundary value problem |?p then v = — log(u) is convex. 

We complete this section by proving Lemma [2j 

Proof of Lemma^ By TheoremJTJ C P (D) = $(it). We integrate by parts and use Q: 

\7u\ 2 dV = [ (Vu,Vu)dV= - [ uAudV 

L> JD JD 

2/p / /■ \ (p-2)/p 

= A / „'•,!]- = A ( J u p ,IV ) { i u p dV 



Rearranging yields 



[ n \Vu\ 2 dV ({ \(p- 2 )/p 
C P {D)= Jd > ' =A( / vPdV) . □ 



ID 



D 
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3 Examples 



We begin with the case of an infinite slab 

S = {(x 1 , ...,!„) : -1 < x n < 1}. 

In order that the variational problem makes sense, one can truncate to obtain S as a limit of Dr 
as R — > oo, where 

D R = {(xi, . . . ,x n ) : -1 < x n < 1, -R < Xj < R}, 

and in the limit we recover the same Euler-Lagrange equation Au + Am p_1 = 0. We look for a 
solution which depends only on x n , which will solve the following boundary value problem for an 
ordinary differential equation (ODE): 



u" + Au p-1 =0, u(-l) = = u(l). 
A quick computation shows 



(14) 



d 
dt 



(V) 2 + |a<i^ = 2u'{u" + AuP- 1 ) = 0, 



so in phase space solutions to ( 14 ) will lie on level sets of the energy function 

E = {u'f + — U p. 

V 



(15) 



Equivalently, for any solution u to ( 14 ) there is a constant E such that 



L 2A 
u' = JE- T vP. 



One can use this last equation to write down all the solutions to (14) up to quadrature, or in 
terms of hypergeometric functions. We sketch some level sets of the energy function ( 15 ) below. 

Indeed, we can compute C P (S) for the slab using just the knowledge that a positive solution 
to ([7]) depends only on x n . On the truncated domain Dr, we have 



f \\7u\ 2 dV = (2R) n - 1 [ (d Xn u) 2 dx n = dST 

JD R J -I 



while 



2/p 



/ r 1 \ 2/p 

vPdV\ = (2i?)i ( ™- 1) [J u p (x n )dx n \ =C 2 R? {n - 1) . 
Taking a ratio we see that $ p (u) = 0(R ( - n ~ 1 ^ 1 ~^), and so 



C P (S) = 



1 <p < 2 

X P = 2 
00 p > 2. 



where in the p — 2 case we have listed the (well-known) value of the first Dirichlet eigenvalue of 
the Laplacian of a slab of width 2. 



Remark 3. Notice that the solution for a slab does not depend on dimension at all. 



G 



10 













I I I I I I I 

i) 5 J 


i i i i 

10 i 


i i i 

15; 


i i i i 


1 1 

20 



Figure 1: Level sets for the energy function (15) with p = 3/2. 



Next we examine the case of the ball. By a moving planes argument IGNNj . any positive 
solution to ([7]) is radial. 

We treat the cases n = 2 and n > 3 separately. First consider the unit disc in the plane. A 
radial solution u will satisfy the ODE 



u" + -u' + KvP- 1 = 0, 



u{l) = 0. 



Change variables by t = — log(r) and let v(t) = u(e ') to get the boundary value problem 

e 2t v + Aw p_1 = 0, «(0)=0, 



(16) 



where the limt-^oo v(t) exists and is a positive number. Here we have used a dot to denote 
differentiation with respect to t and a dash to denote differentiation with respect to r. 
In dimension n > 3 a radial solution will satisfy the ODE 



n- 1 



AuP- 1 = 0, 



u(l) = 0. 



We do a similar change of variables, but this time rescale v by a radial factor, letting u 



-ki 



where again t = — log(r). Again, we will use a dot to denote differentiation with respect to t 
and a dash to denote differentiation with respect to r. Under this change of variables, the ODE 
becomes 



e [2 - k)t {i + (n + 2k- 2)v + k{n + k - 2)v) + Ae 



-k{p-l)t v p-l _ 



0. 



v(0) = 0, (17) 
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where Hindoo e kt v{t) exists and is a positive number. We can eliminate the v term by choosing 
k = and so (17) is now 

2 \ 



e 2 u 



n - 2 \ \ (r,-2)(p-l)t , 

u + Ae 2 v p 1 = o. 



(18) 



In the critical case of p — 5735, the exponential terms coincide, and we obtain the familiar 
differential equation 



n- 2 



v + Aw- 2 = 0, 



which has the energy function 



(19) 



2 2n 

Solutions will lie on level sets of E in the phase plane; see the discussion in [5cj for more details. 



We sketch some level sets of the energy function ( 19 ) below 
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Figure 2: Level sets for the energy function ( 19 1 with n — 3 and p = 6 



4 Comparisons 

In this section we prove some basic comparison principles for minimizers of $ p . The first such 
comparison is domain monotonicity. 
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Proposition 8. If D 1 C D 2 C K ra are bounded domains and p > 1 iftera C p (Di) > C p (Z?2)- 

Proo/. This follows from the fact that L P (D{) n Wl' 2 (D x ) C L P (L> 2 ) n W 1,2 (D 2 ). □ 

Next we fix the domain D and vary p. 

Proposition 9. Fix a smooth, bounded domain D C M. n , and letC p (D) be given by ^ forp > 1. 
TTien i/ie function p i-> C p (D) is continuous. 

Proof. Again, we use the fact that Cq°(D) is dense in L p (D)0 Wq' 2 (D) and take it to be a smooth 
function. In this case, the function 

2/p 



is smooth, so & P (u) is a smooth function of p for fixed « €E C^°(D). The proposition follows from 
the definition of C p . □ 

Our main result in this section is Theorem [3j which reads 

V 2/p C p (D) > V 2/q X q {D) 

for 1 < p < q, where V = Vo\(D). 

Proof. Recall that C^(D) is dense in both L P (D) n W]' 2 {D) and L q (D) n Wq' 2 (Z>), so for the 
purposes of our comparison it will suffice to take u € C^°(D). In particular, u € L P (D) n L q (D). 
Use Holder's inequality on the functions u p and 1, with exponents q/p and — p), to obtain 



2/p / , \2/s 



r — ( I u q dv 



Then, by the variational character of C p and C„, we have 



C, { D) - i„f( ^ |7 ° 1 '":; | > ^ inf( ^'"M . V^CiD), 

which gives the desired inequality. Finally, we can only have equality in Holder's inequality if 
u is constant, which (by the boundary conditions) would force u to be identically zero. This is 
impossible, and the inequality above must be strict. □ 

In dimension two one may take the limit as p — > oo, to obtain 

Coo(^) = inf J £f dV : u E L°°{D) n W 1,2 (D), u £ 1 = lim C P (D). 
{ \\ u \\l™(d) J p ^°° 

By the monotonicty of V 2 ^ P C P (D), this limit exists and is finite. Taking the limit as p — > oo in the 
scaling law (12) with n = 2 shows that CoofV-D) = C oc (D) for each positive r. For a fixed domain 
D, one can find disks rO and KB such that rDCDC KB, so that Coo(i?D) < C^D) < Coo(rD) 
by domain monotonicity. On the other hand, Coo(rO) = Coo(RB) by scale invariance. We see 
that Coo(-D) does not depend on the domain at all, and write Coo for its common value for all 
planar domains. 

We close this section with a preliminary estimate for Coo- 

Proposition 10. Coo < tt- 
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Proof. We may take D to be the disk of radius 1 centered at 0. For any positive 5 < 1, define 
the radial function 

u(r) 

Then ||u||ioom) = 1. Moreover, 

|2, 




D 



1 ' 1 \ 2 , 1 + 6 



so that 



VurdV = 2tt rdr = 7r 



lim 



5-1 J 1-5 



\Vu\ 2 dV 



D 



□ 



5 Extremal domains 

In this section we characterize the domains which are maxima or minima for C p under various 
constraints. We begin with a proof of Theorem |4j that the ball uniquely minimizes C P (D) among 
all domains with a fixed volume. The proof follows the standard proof of the Faber-Krahn 
inequality by symmetrization. 

Proof. Let u be a test function for $ p . Without loss of generality, we can take u £ C^°(D) and 
let 

m = mm\u(x)\. M = max{u(i)j. 

xeD xeD 

For m<t< M let D t = {u > t}. 

Now we define a comparison function : B [m, M] as follows. First let B t be the ball 
centered at the origin with Vol(£> f ) = Vol(Z? t ). Then let u* be the radially symmetric function 
such that B t = {u* > t}. By the co-area formula, 

/ M / ^dr = Vol(A) = Vol(£ t ) = f M f j^-dr. 
Jt JdD T |Vu| J t J 9Bt |Vu*| 



Differentiating with respect to t gives 

dA f dA 



dD t |V«| J dBt |Vu*| 



(20) 



for all t. Then 



f f u p dA , f M „ f dA , 

/) Jm JdD t |VU| J m y ai?t |Vtt| 



as t i v "* i j b 



t p j j^dt = I uldV. (21) 



Now, for m < t < M let 



ijj(t)= \Vu\ 2 dV, ^*(t)= / |Vu.| 2 dV. 

JD t J B t 



By the co-area formula 



ip' = - \Wu\dA, ip'* = - |Vu*|dA 

JdDt JdBt 
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We use the Cauchy-Schwarz inequality, the isoperimetric inequality, and the fact that the normal 
derivative of is constant on dB t to see 

> (Area(<9B t )) 2 = ( / \Vu*\dA 

\JdBt 



dA 



dB t 



We use equation ( 20 1 to cancel the common factor of 

dA f dA 



and so 



dD t |Vu| J dBt |Vu* 



-ip' = / \Vu\dA> / |Vu»|dA = -V'». 

JdD t JdB t 



Integrating this last differential inequality and using ip(M) = = ip*(M) we see that 

\Vu\ 2 dv = ip(o) >v*(o) = / |Vu*| 2 rfy. 



This inequality, combined with (21 1 and ([5]), give the desired inequality on the eigenvalues: 

C P (D) > C P (B). 

Moreover, equality of the eigenvalues forces all level sets dD t to be spheres centered at the origin. 
Also, the equality case of the Cauchy-Schwarz inequality forces |Vit| to be constant on the level 
set dD t . Thus u must be radially symmetric, and so in this case u = u*. □ 

Next we fix the inradius R(D) of the domain rather than the volume, where R(D) is the 
supremum radius of all balls contained in D. 

Lemma 11. Among all bounded domains D C R™ with a fixed inradius, the ball maximizes C p 
for all p > 1 . 

Proof. If the inradius of D is R, then D contains a ball of radius r for each r < R. The result 
now follows from domain monotonicity. □ 

Proposition 12. Let p > 1. Let D be a smooth, convex domain D C M" with inradius R. Let 
u > solve on D and let Um = max{u(x) : x G D}. Then 

2 



9 — n 2AR 

4/< W , (22) 



where 

A.-, * _^ r(1 + '' 



* '„ ^r^- v "r(i + i 



Moreover, (22) is an equality in the case of a slab. 

Remark 4. By the discussion in the introduction, the maximum um is well-defined. 

Hersch [Herl Theoreme 8.1] proved X(D) > tt 2 /(4R 2 ) for the fundamental frequency while 
Sperb |Sp| proved um < R 2 for the maximum value of the torsion function, in each case for a 



convex domain of inradius R. The estimate (22) is a common generalisation of these results. 
More refined results are known in the cases p = 1 and p — 2 (sec MH , for example) . 
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Proof. We follow Section 6.2.2 of [5p]. The P-function 

9A 

v(x) = \Vu(x)\ 2 + —u p (x), 

p 

introduced by Payne, assumes its maximum at the point where u assumes its maximum. Thus 

\Vu{x)\ 2 + — u*(x)<— u p M , 
P P 



which we can rearrange to read 

\Vu(x)\ < 

Let Sd(P) be the distance from the point P where u assumes its maximum to the boundary of 



^yju'u-vrix). (23) 



D and integrate (23) along a line segment which starts at P and terminates at a point on 3D 
closest to P. Then 




The inequality (22 1 follows. Moreover, if D is a slab then (151 implies (23) is actually an equality. 



and so (22) is also an equality. □ 



6 Open questions 

In this final section we collect a small sample of interesting, related questions. 

In the present paper, we have restricted our attention to the minimizer of the the functional 
$ p , which corresponds to the bottom of the spectrum of the eigenvalue equation 

Au + XuP- 1 = 0. 

This functional should have other critical points above the minimizer, which correspond to sign 
changing solutions of the boundary value problem ([7| . What can one say about these higher order 
eigenvalues? If n > 3 and p < (or if n = 2) is the spectrum discrete? Is there a sequence of 
eigenvalues A PJ , with j = 1,2,3,..., such that 

< A P) i < Ap^ 2 < Ap, 3 <■••-> oo? 

We return to the minimum C P (D) = inf{$ p (u)}. In dimension n > 3, does the limit of C P (D) 
exist as p — > from below? It would be nice to characterize the domains on which the limit 
exists, in terms of geometry. We showed in Proposition [9] that the eigenvalue C p is a continuous 
function of p. Is the same true of the eigenfunction? We have shown that, in dimension two, 
lim p _y oo C P (D) = Coo exists and is independent of the domain. Can the bound Coo < tt he 



improved or is it sharp? We proved the inequality (22 ) is realized in the case of a slab. Are there 
any other domains for which ( 22 ) is an equality? 
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